Anisotropic minimal conductivity of graphene bilayers 



OO 
O 

o 

u 

< 



CD 



a 



o 



> 

00 

^' 

O 
00 

O 



X 



Ali G. Moghaddam and Malck Zarcyan 
Institute for Advanced Studies in Basic Sciences, P.O. Box 45195-1159, 45195 Zanjan, Iran 

Fermi line of bilayer graphene at zero energy is transformed into four separated points positioned 
trigonally at the corner of the hexagonal first Brillouin zone. We show that as a result of this trigonal 
splitting the minimal conductivity of an undoped bilayer graphene strip becomes anisotropic with 
respect to the orientation 6 of the connected electrodes and finds a dependence on its length L on the 
characteristic scale I = n/Ak ~ 50nm determined by the inverse of k-space distance of two Dirac 
points. The minimum conductivity increases from a universal isotropic value ct™*" — (8/7r)e /h for 
a short strip L <C £ to a higher anisotropic value for longer strips, which in the limit oi L ^ i varies 
from (7/3)cr""" at 61 = to 3cr^'" over an angle range Ad ~ e/L. 



PACS numbers: 81.05.Uw, 73.23.-b, 73.23.Ad, 73.63.-b 

The recent realization of isolated graphene [l|, l3, Isl , a 
two dimensional hexagonal lattice of carbon atoms, and 
its bilayer [J] has followed by intensive studies which ex- 
plored many intriguing properties of these new carbon- 
based material [5[- They are zero-gap semiconductors 
with their valance and conduction bands touching each 
other at the corners of the hexagonal first Brillouin zone, 
known as Dirac points. This specific band structure 
in connection with pseudo-spin aspect which character- 
izes the relative amplitude of electron wave function on 
two different sublattices of the hexagonal structure, have 
given the carriers a pseudo-relativistic chiral nature [2,lJl- 
The chirality is believed to be origin of most of peculiar- 
ities of quantum transport effects in single and bilayer 
graphene [iiSH- 

One of the most important observations in graphene 
systems is the existence of a nonzero minimal conduc- 
tivity in the limit of vanishing carrier density (at Dirac 
point) [2]. This effect which had been predicted theoret- 
ically long before the experimental synthesis of graphene 
[8|, has been the subject of several recent theoreti- 
cal investigations and in most studies a universal value 
(4/7r)e^//i for minimum conductivity of mono- 



ctq""" = [^/TTje /n tor mmimum conauctivity oi ] 
layer graphene was found 0, [H El, [H, [H Hi. Al- 
though the theoretical value is tt times smaller than the 
value measured in the early experiments, more recent ex- 
periments have confirmed the predicted universal value 
ijmm £pj. .^yj(jg Qj^(\ short graphene strips [l5|, ll6| . 

For bilayer graphene the minimum conductivity is also 
measured to be of order e^/h [J, |l7|. Despite the simi- 
larity of the chiral nature in single and bilayer graphene, 
the low energy spectrum in bilayer is drastically diiTcr- 
ent from the linear dispersion of massless Dirac fermions 
in monolayer. The spectrum in bilayer, which has a 
parabolic form at high energies, acquires strong trigo- 
nal warping at low energies and undergoes the so called 
Lifshitz transition at which the Fermi line is broken into 
four separated pockets [l8|. In the limit of zero Fermi 
energy the pockets shrink into the points of which one 
is located at the Dirac point and three others positioned 
around it in a trigonal form (see Fig. [BJa)). The aim 



of the present letter is to study effect of this Dirac point 
trigonal splitting on the minimum conductivity of the bi- 
layer, which will also allow to distinguish between the 
effects of masslessness and chirality. We employ the real- 
istic model of a wide undoped bilayer strip of length L as 
the scattering region connecting two highly doped regions 
as electrodes. Using a full Hamiltonian which takes into 
account all intra and interlayer hoppings between nearest 
neighbors atomic sites, we find that the effect of the trig- 
onal splitting in the minimum conductivity depends on L 
as compared to a characteristic length £ = n/Ak ~ 50nm 
determined by the inverse of the k-space separation Afc 
of two split Dirac points. 

For a short strip of L '^ £ the effect of trigonal split- 
ting is negligible and cr™™ = (8/7r)e^//i, which is twice 
the minimum conductivity of a monolayer showing that 
in this limit the bilayer behaves as two independent sin- 
gle layer connected in parallel. For finite length strip 
the minimal conductivity increases above cr^™ and finds 
a dependence on the angle 6 between the orientation of 
the electrodes and the hexagonal lattice symmetry axis. 
We find that for a long strip L ^ £ the anisotropic min- 
imum conductivity t7™"'(0) increases from (7/3)(T™™ at 
= to 3tT™™ over an angle range A9 ^ £/L. Our re- 
sults reveals importance of trigonal splitting of the zero 
energy spectrum on the minimal conductivity of bilayer 
graphene. 

To this end, there have been few theoretical investi- 
gations devoted to the minimal conductivity in bilayer 
graphene 0, H HH, [H, 0, ^ . In Ref . ^ a wide bi- 
layer sheet with a constant perpendicular interlayer hop- 
ping is considered to connect two heavily doped electrode 
regions. This model, which ignores the trigonal s plitting , 
results in a minimum conductivity cr™"'- = cr™™ 20|, |2l| , 
establishing on the fact that for high energy electrons 
injected from the metallic electrodes, the constant inter- 
layer hopping does not cause any significant effect and 
the bilayer sheets behaves as two monolayer in parallel. 
On the other hand authors of Refs. [22] and [23( have 
taken into account the effect of strong trigonal warping 
within, respectively. Born approximation and Kubo for- 
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FIG. 1: (a) The bilayer lattice structure projected onto the 
X — y plane. Bonds between Al-Bl (A2-B2) in bottom (top) 
layer are indicated by dashed (solid) lines; Q is the angle be- 
tween the hexagonal lattice symmetry axis and y-axis. (b) 
Constant-energy lines around the corner of the hexagonal first 
Brillouin zone for different energies (measured in units of Lif- 
shitz transition energy e^ ~ ImeV) near the neutrality point. 
Trigonally warped Fermi line is transformed into four sepa- 
rated points with one at the center and three positioned trig- 
onally around it for e = 0. 



mula. They have found an isotropic and constant mini- 
mal conductivity (24/7r)e^//i, which is three times larger 
than a"^'" obtained in Refs. ^ 



2lj . However we note 
that the models employed in Refs. [23, l23[ do not in- 
clude the effect of electrodes which are present in a re- 
alistic experimental setup for conductivity measurement. 
This could be in particular important in graphene con- 
tacts due to the chirality of the carriers and the resulting 
Klein tunneling phenomena p, |7| ■ Our study, taking into 
account both of the trigonal splitting and the electrodes 
effect, reveals anisotropy of the minimal conductivity and 
its dependence on L/£ which also clarifies the origin of 
the disagreement between the two above predictions. 

We consider a ballistic bilayer graphene sheet in x — 
y plane consisting of an undoped strip of length L and 
width W and two heavily doped regions for x < and 
a; > L on top of which bias electrodes are deposited. The 
interfaces between electrode regions and bilayer strip are 
oriented parallel to y axis making an angle 9 with respect 
to the symmetry axis of the bilayer lattice as indicated 
in Fig. [Tfa). 

Bilayer graphene is made of two coupled single layer 
graphene with different sites Al, Bl in the bottom layer 
and A2, B2 in the top layer. The two layers are arranged 
according to Bcrnal stacking in which every Al site of 
bottom layer lies directly below an A2 site in the top 
layer, as shown in Fig. [Ha). Within the tight bind- 
ing model of graphite [5|, [2J] we consider all the nearest 
neighbors hoppings. The unit cell of the bilayer lattice 
structure contains 4 sites Al, Bl, A2, and B2. The intra- 
layer hoppings between the sites Al — Bl and A2 — B2 are 
parameterized by the single energy t w 3eV. There are 
two type of inter-layer hoppings: Al — A2 and Bl — B2 
which are characterized by the energies t± « OAeV and 
is « 0.3ey, respectively. 

The resulting tight binding Hamiltonian is written in 



k-space. The hexagonal first Brillouin zone of bilayer 
graphene contains 6 corners, among them two are in- 
equivalent specifying the valleys K and K'. In the ab- 
sence of inter- valley scattering [g , the valleys are degen- 
erate and it is sufficient to consider only one valley. Low 
energy excitations with 2D wave vector k = (k, q) around 
one of these valleys, say valley K , is governed by Hamil- 
tonian of the form [iSj , 



H(k) 




(1) 



which operates in the space of 4-component spinors of 
the form ^ = 6*'^^+"?'' (-0^11, V'sij '/'A2, "0-82), where each 
component determines the wave function amplitude at 
the corresponding site of the bilayer unit cell. The char- 
acteristic velocities v = 3ta/(2?i),and 113 = 3i3a/(2?i) (a 
is the lattice constant) are associated with the hopping 
energies t and ^3, respectively. We note that the complex 
wave vectors k± = e^*^(fc ± iq) depend on the misorien- 
tation angle 6. 

The quasiparticle spectrum e{k, q) is obtained from the 
eigenvalue equation of Hamiltonian ([1]), which reads 



{e'-\nv\i\'f-e\ii 
= {hv)hj_hv3{e^''^{k- 



2.2 



f inv3kn + \hv3k\H 

iqf + e-'^'^k + iqf). (2) 



It is clear that the right hand side of this equation pro- 
duces trigonal warping of constant-energy lines, with an 
strength given by the ratio 5 = v^/v. Fig. [TJb) shows 
constant-energy lines for three energies e/el = 0, 0.8, 2, 
where El = t^S^ /A is the characteristic energy at which 
the Lifshitz transition takes place. While at e — 2el 
the energy line, despite its strong trigonal warping, is a 
continuous line, at £ = 0.8^^ it breaks into four pock- 
ets whose enclosed area decreases with lowering energy 
and finally at e = shrinks into the four points. The 
central point is located at |k| = and other three leg 
points at a constant distance |k| = Afc = t ±v-i / {Tiv^) 
from the center and in the directions determined by an- 
gles On = arctan(g„/fc„) = 6 + 2{n — l)7r/3, n = 1,2,3 
(see Fig. [lib)). 

Dirac point splitting will have two main effects on the 
transport of carriers through the bilayer strip. First it in- 
troduces a length scale i = tt/A/c = {'K'hv)/{t±5) ~ 50nm 
as the effective range of the scattering potential which 
could mix the states at the four Dirac points. This length 
scale is an order of magnitude larger than interlayer cou- 
pling length l± = hv/t± introduced in Ref. [21|. For 
realistic graphene samples of few 100 nm length, L > i 
and the variations over this new length scale should be 
taken into account. The potential profile varies over the 
length L of the strip. For L '^ £ the scattering potential 
is short enough to cause strong inter Dirac point scatter- 
ing and consequently the scattering states at these points 



are mixed. On the other hand for L ^ £ the states of 
the Dirac points are well separated and do not mix. Sec- 
ondly, it causes an anisotropy of the scattering states due 
to the orientation of the leg points. As we will explain in 
the following, these effects result in a length-dependent 
anisotropic minimal conductivity for the bilayer strip. 

Within the scattering formalism we find the transmis- 
sion amplitude of electrons through the bilayer strip. An 
electronic state is specified by the energy e and the trans- 
verse wave vector g, which are conserved in the scattering 
process. We find the eigenstates of Hamiltonian ([T|) in 
three regions of left (a; < 0) and right {x > L) electrodes 
and the bilayer strip (0 < x < L). In general for a given 
e and q there are 4 values of longitudinal momentum k 
(solutions of Eq. ^) with 4 corresponding eigenstates 
in each region. 

Inside the strip at the neutrality point e = 0, the 4 
solutions (fci, z = 1, ..., 4) of Eq. ([2]) have the form 



fci,: 



Ka = iq + 



\i±\i + i 



o3ie 



2e ■ \ TT 

The corresponding eigenstates are given by 



)• (3) 



03,4 



Jki^-iX+iqy 



Jks^iX+iqy 



(0,- 



l± 



-,ki^2-iq,0), (4) 



(fc3,4 + ig,0,0,-T— ). 



(5) 



In general for a given q all 4 states inside the strip 
are evanescent having complex ki{i = 1,...,4) with ex- 
ceptions of the Dirac points with qo = and g„ = 
(tt/^) sin(6',i) (n = 1,2,3) at which two of kiS are real 
representing propagating states in the strip. 

Inside the highly doped electrode regions a very large 
potential —Uq is applied. We assume that for all states 
contributing in transport hvq <i; Uq provided that Uq be 
much larger than all other energy scales in Hamiltonian 
([T|) . By this assumption the longitudinal momenta of all 
states inside electrode regions have a constant magnitude 
^0 = Uq/Kv ^ q. Inside each electrode for certain q there 
arc two right-going eigenstates. 



lR _ pikox+iqy ( 



.±-c -'(l,e-'^±l,±e^''), (6) 

and two left-going eigenstates. 



e-*'^+*«^(l, ~e"\ ±1, Te''"). 



(7) 



For two left-going incident states from the left elec- 
trode {x < 0) the scattering states in three different re- 
gions have the form 



A^ _1_ -r-^J^L _i_ ^±^L .T < 



*^ 



where the coefficients cf and reflection and transmission 



< a; < L, (8) 

X > L. 



amplitudes ri , ti have to be determined by imposing the 




FIG. 2: Minimal conductivity in units of cr^^" = (8/7r)(e^//i) 
of a wide strip of undoped bilayer graphene versus its length 
L for different orientations 6 of the hexagonal lattice sym- 
metry axis with respect to the electrodes. L is measured in 
units of ^ = Tv/Ak with Afc being the k-space distance of two 
trigonally split Dirac points. The isotropic minimal conduc- 
tivity a^^" at L <^ £ increases with L/i to higher values and 
becomes anisotropic. 



continuity condition of the wave functions at the bound- 
aries X ~ 0, L. 

We calculate the conductance at zero temperature 
from Landauer-Buttiker formula, 



G W f°° 
TT^TT T{q)dq, 



(9) 



where T{q) = |4P + \t+? + 1*1 P + \t-? is the sum 
of transmission probabilities of the two states ^±, and 
Go = Ae^ /h is 4 times of conductance quantum to take 
into account the valley and spin degeneracies. The con- 
ductivity of the bilayer strip is obtained by the relation 
(J = {W/L)G. 

We have obtained T{q) as a function of the length L 
and the orientation angle 0. For a short strip with i <C ^ 
the transmission probability takes the form 



T{q) 



1 



1 



cosh^ [{q — qc)L] cosh^ [{q + qc)Ly 



(10) 



which shows two maxima at the point q = ±qc = 
±arcsinh(L/2/j_)/L. Around these points T(g) decays 
exponentially within a scale of order Aq ~ 1/L which for 
a short strip is much larger than Afc = ir/i. Thus T{q) is 
almost constant within the scale Afc which implies that 
the Dirac point splitting and the associated anisotropy 
of the spectrum are not revealed in the transmission of 
the carriers. This is the result of strong mixing of the 
states around the 4 Diract points via scattering through 
the bilayer strip. The resulting isotropic minimum con- 
ductivity a^™ = 2(70"*" which is twice of the minimum 
conductivity for a single layer. 

For a finite L/i the minimum conductivity increases 
above cr™*" . This is shown in Fig. [5] where we have 



plotted a'^'^^ as a function of L/£ for different orientations 
6. For L/i ^ 1 the increased minimum conductivity finds 
a 6'-dependence as the result of trigonal splitting of the 
Dirac point. In this case the anisotropy of a"^'-^ is spread 
over the range < 9 < tt/6. We note that the increase 
in cr™™ is associated with an oscillatory variation due to 
quantum interference effects in bilayer strip. Increasing 
L/£ further to approach the limit of a long strip L ^ £ 
the range of anisotropy becomes narrower. In this limit 
and for the angles not too close to 9 = 0, 7r/6 we obtain 
the following relation for the transmission probability 



n<i) = E 



1 



„f^3COsh2[a„(^)(g-q„(0))L]' 



(11) 



where the summation is taken over 4 transverse co- 
ordinates of the Dirac points, Qn, with ap = 1 and 
ani9) = 3/(5 -4cos[20„]) for n = 1,2,3. This result 
shows that the transmission probability consists of 4 res- 
onant peaks at the points g„, whose width is of order 
Aq ^ 1/L << Afc. The effect of the Dirac point split- 
ting is, thus, revealed in the transmission process. For 
9 approaching and tt/6 the transverse coordinates of 
the two resonant peaks qq -^ qi and qi -^ (72, respec- 
tively, and the corresponding peaks overlap. For these 
cases Eq. (fTTj) is not applicable, since we have consid- 
ered four Dirac points contribution independently in this 
equation. In the case oi 9 = from Eqs. ^ and ^ 
we find that the contribution of the conductivity from 
the peaks at go = 91 = is of order S^e'^/h which is 
negligibly small. So there are only contributions from 
the points <3'2,93 which results in a minimum conduc- 
tivity CT^!!JJ = (7/3)cr^"\ In contrast to 6* = case, 
for 9 — t: /Q the overlapped resonant peaks makes the 
same contributions as two independent peaks. This gives 
a minimum conductivity a™'™ ~ 3cr™"\ Wc find that 
for a long strip this value of the minimum conductiv- 
ity is valid for all orientations ^ 6* < 7r/6 except of 
the angles 9 < £/L where the two peaks goi^i have a 
significant overlap. Thus the minimum conductivity is 
anisotropic over the range /S.9 ^ £/ L and an amplitude 
Act™*" = (2/3)ctX"". 

In conclusion, we have studied conductivity of a wide 
strip of undoped bilayer graphcne which connects two 
highly doped electrode regions. We have shown that due 
to the trigonal splitting of Dirac point at zero Fermi en- 
ergy, the minimal conductivity a™"™ of the strip finds a 
dependence on the lattice symmetry axis orientation 9 
with respect to the electrodes. The anisotropy of cr""" 



depends on the length of strip L as compared to the char- 
acteristic length scale £ = ir/Ak ~ 50nm determined by 
the inverse of the k-space separation of two Dirac points. 
For a short strip of L ^ ^, a"^'^'^ takes an isotropic uni- 
versal value a^™ = (8/7r)e^//i. For longer strips the 
minimal conductivity increases above this value in an 
anisotropic way. We have found that in the limit of 
L ^ £ the anisotropic minimal conductivity grows from 
(7/3)cr^™ at 6* = to 3cr^" when orientation is changed 
by the angle A9 ^ £/L. 
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